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1. Introduction 

The ^^C(a,7)'^0 reaction has been considered to be one of the most important nuclear reac¬ 
tions in the nucleosynthesis of elements in the universe [1]. However, due to the Coulomb barrier, 
the cross section at the center-of-mass energy Ec,m. = 0.3 MeV corresponding to helium burning 
temperatures is too small to measure using present laboratory technologies. In this circumstance, 
the low-energy cross section is extrapolated from the available experimental data with theoretical 
guides, and it is converted into the nuclear reaction rates for astrophysical applications. 

The representative studies of the ^^C(a,7)^^0 reaction rates have been provided by KU02 
[2] and BU96 [3], based on the /^-matrix analyses. For the extrapolated astrophysical S'-factor at 
Ec.m. = 0.3 MeV, which is used instead of the cross section, KU02 estimates Se\ = 76 ±20 keV b, 
Se 2 = 85±30keVb, and Scasc = 4±4keVb for theSl, £"2, and cascade transitions. BU96 predicts 
‘^£1 = 79 ± 21 keV b, Sei = 70 ± 70 keV b, and S'casc = 16 ± 16 keV b. The derived reaction rates 
have been reported to give a different temperature dependence from that in the pioneering Fowler 
compilation of the reaction rates (CF 88 [4]). The so-called NACRE compilation [5] following 
CF 88 is the milestone of the compilation work, providing the evaluated rates to the astrophysical 
community. 

In the present study, we illustrate the calculated results of the i'-factor and 7 -ray angular dis¬ 
tribution, and we compare the derived reaction rates (KA12) [ 6 ] with the previous studies (CF 88 , 
NACRE, KU02, BU96). We use the potential model [ 6 , 7]. This is the one of direct nuclear reac¬ 
tion models, which describe the reaction process, in which a few degrees of freedom of motion are 
activated. Einally, we show an analytic expression of our reaction rates. 

2. Potential models 

We first describe a brief outline of our direct capture potential model [ 6 , 7]. The radiative cap¬ 
ture cross section a is calculated from a oc \ {(pf\Mf_\(pi)\'^, where 9,, (pf, andM^ are the scattering 
wave, the bound state, and electromagnetic operator, respectively. A denotes multipolarity of the 
transition. The input parameters of the model come from the potential used to generate 9 , and 
(Pf. The effective chai'ge in Mf is also used as an adjustable parameter to minimize the difference 
between the theoretical values and the experimental cross section data. The astrophysical i'-factor 
is defined as Sex = ^c.m. exp(27rt])cJ, where 77 = I2e^/hv. v is the relative velocity. 

The (pi is generated from the potential that describes phase shifts and differential cross sec¬ 
tions for a-i-'^C elastic scattering [ 6 , 7, 8 , 9]. The internuclear potential should be chosen in the 
appropriate strength. To obtain the correct strength, we refer to nuclear rainbow scattering at high 
energies [10, 11, 12]. The potential used here resembles the double-folding potential with M3Y 
interaction [13], and it produces the a+^^C elastic cross sections shown in Eig. 1. Eor the bound 
states (pf, we classify the states into two groups. One is the states well described by a-cluster 
models. (Oj: £x = 6.05 MeV, and 2^: Ex = 6.92 MeV.) E^ denotes the excitation energy of *^0. To 
make the a-cluster states, we use the potential for the a-i-^^C rotational bands [7]. Another is the 
shell model states. (Oj^: ground state, 3]^: Ex = 6.13 MeV, and 1]“: Ex = 7.12 MeV.) To make these, 
we use the so-called separation energy method, in which the potential is adjusted so as to repro¬ 
duce the experimental separation energy. This is the same as a textbook method in direct reaction 
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Figure 1: The excitation functions of a+*^C elastic scattering below Ec.m = 5 MeV [7]. The solid curves 
are the calculated results with potential scattering. The arrows indicate the resonance energies of the states 
in the rotational bands. The experimental data are taken from [8]. 


models, e.g. distorted wave Born approximation (DWBA) [12]. It gives appropriate wavefunctions 
in the peripheral region that nuclear reactions are sensitive to. The large difference of the potential 
between the bound state and scattering state may indicate that the wavefunction is not appropriate 
in the assumed configuration. However, it is good for a description of direct reactions as a doorway 
to making fused nuclei from scattering states. 

The Maxwellian-averaged reaction rates Na (cjv) [1, 5] are calculated from the theoretical cross 
section by the potential model. 

3. Results 

Before comparing the reaction rates with CF88, NACRE, KU02, and BU96, we show the 
theoretical calculation of the S'-factor and /-ray angular distribution from the potential model. 

The calculated S'-factor for El and E2 transitions of the ^^C(a,7o)'^0 reaction is compared 
with the experimental data [2, 14] in Fig. 2. From this figure, the E2 transition is found to be 
enhanced at low energies. The /-ray angular distribution of the ^^C(a,/o)'^0 reaction is shown in 
Fig. 3. The potential model (solid curves) appears to reproduce the trend of the peak and valley 
of the /-ray angular distribution. At Ec.m. = 2.267 MeV, the /-ray angular distribution has a single 
peak, that means the p-wave dominates the transition. At low energies, the d-wave seems to be¬ 
come more important. Although the E 1 component appears to deviate from the experimental data 
in Fig. 2(a), the angular distribution seems to be reproduced by the potential model. The uncer¬ 
tainties in the experimental /-ray angular distribution have recently been discussed in [15]. The 
extrapolated values of the S'-factor at Ec.m. = 0.3 MeV are listed in Table 1. The total S'-factor from 
the present work is consistent with KU02, NACRE, BU96. In contrast, our result for the E 1 and E2 
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Figure 2: The astrophysical 5-factor of the 

^^C(a,7o)'®0 reaction; (a) £1 and (b) E2. The 
solid curves are the results obtained from the po¬ 
tential model. The experimental data are taken from 
[2, 14], 



Figure 3: The y-ray angular distribution of the 
*^C(a,7o)*®0 reaction [6]. The experimental data 
are taken from [2, 14]. 


Table 1: Astrophysical 5-factor at = 0.3 MeV. The 5-factor is listed in keV b unit. The total 5-factor 
includes the contribution from the cascade transition through the excited states of ^®0. 



This work 

KU02 [2] 

NACRE [5] 

BU96 [3] 

El 

E2 

Cascade 

Total 

3 

150+|‘ 

18±4.5 

171-t46 

76 ±20 

85 ±30 

4±4 

165 ±54 

79±21 

120 ±60 

199 ±81 

79±21 

70 ±70 

16±16 

165 ± 107 


components is different from the previous studies. However, let us recall that the decomposition of 
El and E2 is performed from the analysis of the y-ray angular distribution. 

The subthreshold 2^ state and the resonant 1 2 state at E^ ^9.6 MeV are the members of the 
a-i-^^C rotational bands in [6, 7]. They appear to have the broad width of electric transition, as 
shown in Fig. 2. The narrow resonances do not interfere with the present result so much, and they 
may be negligible [16]. The /^-matrix may not give a good description of the *^C(a,7o)'^0 reaction 
at the astrophysical energies, because it cannot explain e.g. the weak coupling of the system, the 
state with a broad width, and large violation of the isospin selection rule [16]. The experimental 
studies of the total cross sections [17] have also been performed, recently. 

The asymptotic normalization constant (ANC) can be obtained phenomenologically, = 
5.00 X 10^^ fm^* for Ij", = 2.03 x 10^° fm^' for 2^. These values are consistent with the 
previous studies [18, 19]. 

In Fig. 4 the derived reaction rates (KA12) are compared with those from the representative 
studies. The panel (a), (b), (c), and (d) show the ratio of the reaction rates to CF88, NACRE, 
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Figure 4: Comparison between our reaction rate KA12 and one from the previous studies. The reaction 
rates are expressed in the ratio to (a) CF88, (b) NACRE, (c) KU02, and (d) BU96. The solid curves are 
the recommended reaction rates. The shades are uncertainties estimated from the model parameters. The 
dense and thin shades come from the variation of the potential parameters and the strength of the cascade 
transition. The dotted and dashed curves are the upper and lower limits estimated by the respective group. 


KU02, and BU96, respectively. We find from this figure that KA12 seems to be close to BU96. 
The theoretical uncertainties of the reaction rates are estimated from the variation of the model 
parameters [6]. The dark area comes from the variation of the potential parameters; the thin hatch 
comes from the strength of the cascade transition of ^^0. 

Our reaction rates below Tg = 3 (KA12) are approximately given in the simple expression in 
unit of cm^mol^^s^^ 


Na{ov) = 1.16xl0^rg ^'^^exp 


32.369 

rpl/?! 


T9 

3.17 


1 - 0.371 Tg - 0.106 + 0.26479^, 


where Tg is temperature in units of 10^ K. The upper and lower limits, Na{ov)u and 77 a(cjv)l, are 
expressed by 

2 " 


4/3 

exp 

32.330 

r 79 


^2.80 

4/3 

exp 

32.363 

( 79 

y,l/3 

^4.04 


1 - 0.0686 Tg - 0.845 Tg^ + 0.711 Tg^ , 


1 - 0.422 Tg + 0.124 + 0.0894 . 


The tabular form of the reaction rates can be found in [6]. 


4. Summary 

The astrophysical S'-factor of ^^C(a,7)^^0 has been investigated with the potential model. The 
derived reaction rates have been compared with those from CF88, NACRE, KU02, and BU96. 

The total S'-factor obtained from the potential model is concordant with NACRE, KU02, and 
BU96. However, the E2/E1 ratio is different from that of the previous works. Erom the potential 
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model, the S'-faetor at low energies is predieted to be dominated by E2 transition to the ground 
state. The E\ and E2 S'-faetors at Ec,m. = 0.3 MeV are Sei ~ 3 keV b and Sei = 150ljy keV b, 
respeetively. The sum of the easeade transition through the exeited state of is ^casc = 18 ±4.5 
keV b. The y-ray angular distribution seems to be reprodueed with the potential model, whereas 
the E\ S'-faetor at low energies appears to deviate from the experimental value. We do not invoke 
a eompensatory large eontribution beeause of the weak eoupling feature of the system. 

The derived reaetion rates at low temperatures seem to be eoneordant with CF88, NACRE, 
KU02, and BU96. For astrophysieal applieations, our reaetion rates below Tg = 3 have been pro¬ 
vided in the analytie expression. 
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